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Abstract 

The theoretical description of compact structures that share some features with mass varying 
particles allows for a simple analysis of the equilibrium and stability for massive stellar bodies. We 
investigate static, spherically symmetric solutions of Einstein equations for a system composed by 
nonbaryonic matter (neutrinos or dark matter) which forms stable structures through attractive 
forces mediated by a background scalar-field (dark energy). Assuming that the dark matter, or 
massive neutrinos, consist of a gas of weakly interacting particles, the coupling with the scalar 
field is translated into an effective dependence of the mass of the compounding particle on the 
radial coordinate of the curved spacetime. The stability analysis reveals that these static solutions 
become dynamically unstable for different Buchdahl limits of the ratio between the total mass- 
energy and the stellar radius, M/R. We also find regular solutions that for an external observer 
resemble Schwarzschild black-holes. Our analysis leaves unanswered the question whether such 
solutions, which are both regular and stable, do exist. 
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I. INTRODUCTION 



The search for the mechanism responsible for the onset of the accelerating phase of the 
Universe continues to stimulate interesting and fascinating discussions on many key issues in 



astrophysics and cosmology 



In the context of models account for the accelerated 
expansion of the Universe, a nontrivial coupling between the nonbaryonic matter sector and 
the sector responsible for the acceleration of the Universe has often been considered 

aaa. 

A possible trigger for the acceleration arises through the interaction of the quintessence field 



identified with dark matter 



llj . neutrinos 
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or is unified 



with dark energy as in the Chaplygin gas model 

Generically speaking, nonbaryonic matter interacting adiabatically with a background 
scalar field gives origin to the so-called mass varying mechanism, originally conceived to 



address the dark matter issue, and later reformulated to study cosmological neutrinos 
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171]. In its simplest realization, it consists in considering a scalar field effective potential 



with an implicit dependence on the neutrino mass. Due to the constraint imposed by the 
scalar field, the mass eventually grows close to its present value, when neutrinos form a 
nonrelativistic (NR) fluid and the interaction with the scalar field locks the mass evolution. 
The potential energy of the dark energy component becomes then the dominant contribution 
to the Universe's energy density and the cosmic acceleration ensues. However, in models 
with increasing particle mass, the homogeneous cosmological solution is usually unstable to 
perturbations. In fact, several forms of coupling of mass varying particles to quintessence 
fields lead to instabilities associated with an imaginary speed of sound for the fluid in the NR 
regime. Such instabilities result in the exponential growth of cosmological perturbations 25|. 



The Universe becomes then inhomogeneous and these particle lumps form denser structures. 

In this work we quantify the equilibrium and stability conditions for static, spherically 
symmetric objects that share some features with cosmological scenarios with mass varying 
particles. These exotic astrophysical objects are compact lumps of nonbaryonic matter, held 
together by gravity and the attractive force mediated by the background scalar field. The 
static mass varying behavior is achieved from the dependence on the scalar field radial co- 
ordinate. In our analysis, we shall consider that such stellar structures are described by the 
Schwarzschild's solution of a sphere with constant energy density, p, and pressure, p, which 



2 



drops from its central value to vanish at the boundary. We shall determine the equilibrium 
conditions between the total mass-energy, M, and the spherical radius, R. Actually, in our 
approach the variable mass behaves like an additional extensive thermodynamic degree of 
freedom, due to its dependence on the scalar field, being therefore like a chemical potential 
from the thermodynamic point of view. Under adiabatic conditions, the thermodynamic 
pressure and its explicit dependence on r are given in terms of p{r), and can be explic- 



itly computed 
equilibrium 
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rom the Tolman-Oppenheimer-Volkoff (TOV) equations for the hydrostatic 



271]. This approach justifies the assumption of a specific analytical depen- 



dence of the variable mass on the radial coordinate of the massive object, m{r), at a given 
distance from the center. The stability analysis reveals that our static solutions become dy- 
namically unstable for different Buchdahl limits of the ratio between total mass-energy and 
stellar radius, M/R. Assuming that such compact structures may describe massive neutrino 
lumps, we also verify how the structure formation could eventually modify the cosmological 
predictions for the neutrino masses. In particular, we find that there are some physically 
acceptable solutions which give rise to stable lumps. Our analysis leaves open the issue 
related with stable regular solutions that closely reproduce the conditions for the formation 
of Schwarzschild black holes. 

This work is organized as follows. In section II, we review well-known thermodynamic 
properties for neutral matter in order to properly formulate a form of the first law of ther- 
modynamics that is consistent with the mass varying mechanism. In section III, we review 
the fundamental aspects of the mass varying mechanism in the Friedman-Robertson- Walker 
(FRW) cosmological scenario. In the context of structure formation, we establish the con- 
nection among the energy density p, the variable mass m and the particle density n, for 
an adiabatic system of mass varying particles in a compact structure. The corresponding 
equilibrium conditions and the mass defect are determined in section IV. In section V, we 
perform the analysis of stability of the equilibrium configurations and obtain the conditions 
for neutrino lumping. We also verify how the cosmological predictions for neutrino masses 
are modified. Finally, we summarize our findings and discuss their implications in section 
VI. 
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II. GENERAL THERMODYNAMIC PROPERTIES FOR NEUTRAL MATTER 



In order to analyze the cosmological problem and the theory of stellar evolution, it is 
necessary to get information on the interactions among the particles which make up the 
astrophysical bodies. Regardless the specific nature of these forces, a common property is 
the additivity of the interaction energy for a macroscopic system: if the system is divided into 
macroscopic parts, then the interaction energy between these parts will be negligibly small, 
and for short-range forces we can introduce the concept of specific energy e = E/N = p/n, 
where p is the energy density, the total energy, E, in a volume, V, and n is the density of 
particles to which we relate the energy and is the total number of particles, = nV. 

Furthermore, the main manifestation of short-range interactions in a macroscopic system 
is a nonvanishing pressure p. The pressure is the quantity which allows for defining the 
interaction force between two parts of a system and depends only on the state of the matter, 
p = p{n, s), i. e. on the density of particles n and on the specific entropy s = S/N. Finally, 
we have to introduce other two thermodynamic intensive variables: the temperature, T{n, s), 
and the chemical potential, p{n, s). 

The macroscopic version of the first law of thermodynamics leads to the familiar expres- 
sion which relates the abovementioned quantities 

d^ = -p{V, S) dV + T{V, S)dS + J2 H dA^,', (1) 

i 

where the indices j correspond to different types of particles. Equally fundamental is the 
conservation of total number of particles. Suppose a fluid element V whose moving walls 
are attached to the fluid so that no particles flow in or out. As the fluid element moves, its 
volume changes, but the total number of enclosed particles remains fixed, 

d(E^^) = 0' (2) 

i 

which can be expressed as 

dN = d{nV) =0, (3) 

for neutral matter, i. e. for fiuid made up by a single type of particle. In this case, the 
volume occupied by a particle can be expressed as = 1/n, the corresponding entropy by 
S = N s = nV s = and the thermodynamic relation from Eq. ([T]) reads, 

de = d = -p{n, s) d + T{n, s) ds, (4) 



which sets 



's^j ^^i^^j dn + nT{n, s)ds. (5) 



?! 



The analysis further simphfies the fluid, already assumed to have the same composition, has 
the same entropy per particle, s, and is in a state of adiabatic flow, so that pressure and 
energy density are related by 

(^) = /'('^' ■^)- (^) 
The influence of the interactions on the thermodynamics and the assumptions about 
the spectrum of rest masses in the context of neutral matter stellar structures have been 
recurrently examined the literature. In the following analysis, we shall consider neutral 
matter with no color, electroweak and electric charges such as (cold) dark matter and NR 
neutrinos. These considerations will allow us to describe interactions of neutral particles 
within stellar objects whose evolution is driven by the coupling to the dark energy scalar 
field. 



III. THE MASS VARYING MECHANISM 



In the so-called mass varying mechanism 12|, ll3|, [ij, l28| a coupling between a relic particle 



and a light scalar field is introduced. This fluid is identified with dark energy or dark energy 
plus dark matter [l^, [l^ . As a consequence of this new interaction, the mass of the relic 
particle is generated from the vacuum expectation value of the scalar field and becomes 
linked to its dynamics, that is m = m(0) for such a relic particle. 

In the usual MaVaN framework, neutrinos are coupled to a light scalar field which is 
identified with the dark sector. Presumably, the neutrino mass rriiy has its origin on the 
vacuum expectation value (VEV) of the scalar field and, naturally, its behavior is governed by 
the dependence of the scalar field on the scale factor. Given a particle statistical distribution 
f{q), where g = Tg being the neutrino background temperature at present, in the flat 
FRW cosmological scenario, the corresponding energy density and pressure can be expressed 
by 

+ (7) 



where sub-index denotes present-day values, such that = 1. It can be easily shown that 

dp{a, 0) 



m{(t))- 



{p{a,cj)) - 3p{a,(f>)). 



(8) 



dm{(t)) 

From the dependence of p on a, one can obtain the energy-momentum conservation for the 
mass varying fluid, 

dm{(j}) dp{a, 4)) 



p{a, <j>) + 3H{p{a, (/)) + p{a, (/))) = 0- 



(9) 



dcp dm{4>) 

where H = a/ a is the expansion rate of the Universe and the dot denotes differentiation 
with respect to cosmic time (' = d/dt). It is important to emphasize that the coupling 
between relic particles and the scalar field as described by Eq. (IHl) is restricted to times 
when they are NR, i. e. ^'^^"'/f) — n{a) oc [12I 



Ij, l28|. In opposition, as long as particles 



are relativistic (T(a) = T^o/a » m{(t){a)))^ the decoupled fluids should evolve adiabatically 
since the strength of the coupling is suppressed by the relativistic increase of the pressure 
(p ~ 3p). The mass varying mechanism is essentially driven by Eq. (jH]), which translates the 
dependence of m on into a dynamical behavior. 

However, the mass varying mechanism can play a salient role in the structure formation. 
It can be assumed that cosmological neutrinos, dark matter, or some generic form of non- 
baryonic matter, compose a gas of weakly interacting particles. The mass of these particles 
depends on the value of a slowly varying classical scalar field. This dynamics is similar to 



that of the so-called cosmon fields which can lead to the formation of cosmon lumps 
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29| 



Cosmon lumps are bound objects, for which the scalar and gravitational fields combine 



to form nonlinear so 
detail in Refs. 30 



utions of their field equations. Such solutions have been discussed in 



31 



32( 1 for situations where the infiuence of the cosmon potential is 



negligible. Their characteristic features are: (i) lumps are designated by their mass and a 
suitably defined cosmon charge; (ii) the black-hole event horizon disappears for all solutions 
with nonzero cosmon charge; (iii) for small cosmon charge, the solutions closely resemble 
black-hole solutions when seen from the outside, and the horizon is replaced by a sharp 
transition to a region of high red-shift. 

For classical particles, the action of the cosmon coupled to mass varying particles {v or 
DM) can be written as 31| 



S 



+ \g^''d^(l)d,(l) + ¥{<!>) + p{ct>) 



(10) 



where G is the Newton constant, R is the Ricci curvature and p{(j)) is the corresponding 
mass varying particle density. For static, spherically symmetric solutions, one employs 
the Schwarzschild metric given by ds^ = —B{r)dt'^ + A{r)dr^ + (d^^ + sm^{6)d(f^), where 
A = B^^ = (1 — 2G'M/r)^\ for which we set G = 1, and r is the radial coordinate. From 
the field equations we can write 



dr" 



2 B' 



2A) dr 



^ f dV{4>) dp{4>) \ 
V d0 d(p ) 



from which the explicit dependence of on r is obtained. 



Assuming the adiabatic approximation at cosmological scales [25|], the cosmological sta- 



tionary condition 



13|, 



14 



16| can be applied to Eq. ffTTl) in order to suppress its right 



hand side. The analytical dependence of the scalar field on the radial coordinate effectively 
disappears for large values. Faraway from the concentration of matter, i. e. for r larger 
than the compact structure radius, one recovers the background of a flat FRW Universe. 
In this case, considering back the usual cosmological time-dependent terms, the Eq. (fTTj) 
should reduce to the energy conservation equation. 



, orri , , dlnm(0) ^ 

+ '^'^+i^ + ^^^^-'^^ 



0, 



(12) 



and the stationary condition could be equivalently assumed for the static approach, i. e 

dV{(f)) dlnm(0) 



+ 



-(p-3p) = 0. 



(13) 



This equation does not depend on the relation of (f) with spacetime coordinates or metric 
components, that is, the radial coordinate, r, of a curved region or the scale factor, a. For this 
reason, the simplest scenario is obtained from the extension of the cosmological stationary 
condition to a static conflguration in Eq. (ITTl) . The adiabatic regime depicted from the 
gradual change expressed by Eq. (fT3ll can be compared to the minimum of the effective 



potential of a chameleon models 



34j . where the chameleon scalar fleld are modifled 



through spatial variations in the interior of massive bodies. In addition, the corresponding 
energy momentum conservation law is exactly what one encounters in the Einstein frame of 
an scalar-tensor gravity model 35|, ISGj. In Ref. 251], this adiabatic regime and its domain 
of validity is reviewed in a general context. 
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Introducing the abovementioned Schwarzschild expressions for A{r) and B{r), and assum- 
ing the stationary condition, Eq. (|TT|) can be simphfied as 



dr^ 



+ 



r -2{M/R){r^/R^)\ dr 
d^ 



^^0 



dr2 



+ 



2(r - M)\d(j) _ ^ 



which yields the solution 



° 2M 



In 1 



2M 



0(0 



r - 2M y dr 

(r < i?), 
1 



(r < R), 
(r > i?), 



(14) 



v/(2M7iF 



arc 



tanh(v/(2M/i?3 



(r > i?),(15) 



where 0o'i""'' are constants to be adjusted in order to match the boundary conditions. In 
the Newtonian limit where A{r) ^ B{r) ^ 1, the constants match one each other so that 
()o,i and the above solutions reduce to 

01 



0,1 



Lout 
^0,1 



(r) = 00 + 



(16) 



Given that depends on the radial coordinate for generic curved spaces, we quantify the 
changes due to our approach through some mass dependencies on r as illustrated in Fig. [H 
Clearly, the form of rri(0) and of the relation between pressure, p{r), and density, p{r), may 
lead to quite different scenarios. 

In the pictorial representation of our results we denote y = r/R, where R is the radius of 
the spherical lump. Any prescription for m{r) is model dependent, i. e. arbitrary functions 
for 'm{r) are equivalent to arbitrary functions for 0(r) and m{(f>). In particular, the stationary 
condition approximation justifies the analytical dependence given by Eq. (fT8|) . However, 
other approximations could be adopted. 

Let us now turn to scenarios involving neutrinos. Despite its impressive phenomenological 
success, it is widely believed that the Standard Model (SM) of particle physics is actually 
only a low-energy approximation of an underlying more fundamental theory. In this respect, 
the interplay with the cosmology can be an important guideline to obtain insights on the 
nature of the more fundamental theory. Turning to the subject of neutrinos in the SM, the 
most natural way to explain the smallness of their masses is through the seesaw mechanism, 
according to which, the tiny masses, m^, of the usual left-handed neutrinos are obtained via 
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a very massive, M , sterile right-handed neutrino. The Lagrangian density that describes 
the simplest version of the seesaw mechanism through the Yukawa couphng between a hght 
scalar field and a single neutrino fiavour is given by 



^ = rriLRi^L^ji + M{(j))uRijR + h.c. 



(17) 



where it is shown that at scales well below the right-handed neutrino mass, the following 
effective Lagrangian density emerges [371] 



m 



LR 



M{<l>) 



J^lJ^r + h.c. 



(18) 



Phenomenological consistency with the SM implies that logarithm corrections to the above 
terms are small. Furthermore from solar, atmospheric, reactor and accelerator neutrino os- 
cillation experiments it is gathered that neutrino masses given by m,^ ~ m|^/M(0) lie in the 
sub-e\^ range. It is also clear that promoting the scalar field (p into a dynamical quantity 
leads to a mechanism in the context of which neutrino masses are time-dependent. Asso- 
ciating the scalar field to the dark energy field allows linking NR neutrino energy densities 



to late cosmological times 



3, 13, Q: 



15 



161] . This scenario was implemented through the 



adiabatic (stationary) approach via Eqs. (fTTjl . ( fT2i) and ( fT3|) . It is evident that this approach 
is fairly general, and that the same analysis can be performed in the framework of mass vary- 
ing particles coupled with cosmon fields in curved spacetime. Analogously, an approximate 
solution for a spherically symmetric compact object of radius R with homogeneous density 
were already obtained in the quoted chameleon framework 33] 



A. Thermodynamics of "noninteracting" mass varying particles 

To describe the connection among the extensive quantities p, m and n in an adiabatic 
system of mass varying particles in a star (or any stellar object), the relevant thermodynamic 
equations are summarized by Eqs. (I6l)-(IH]). In fact, in order to get a closed system of 
equations, an equation of state that relates the pressure in terms of the energy density 
and the specific entropy, p = pip,s). Moreover, due to the adiabatic conditions set up, 
the pressure and its explicit dependence on the radial coordinate, r, is simply given in 
terms of p(r), which can be explicitly obtained from the TOV equations for the hydrostatic 
equilibrium. 
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By assuming that the hydrostatic pressure, p = p{p{r)), leads to a univoquous dependence 
of p on the radial coordinate for a symmetrically spherical distribution of matter, we have 
to obtain a density p{m{r),n{r),r) which satisfies the system of partial differential equations 
given by 

dp dp dm ^ dp dn ^ dp 
dr dm dr dn dr dr ' 
and Eqs. dHD-dHD- Eliminating p from Eqs. (JH])-®, we obtain 

4p = m|^ + 3n|^, (20) 
dm dn 

and consequently, 

p(m,n) = Km^-'"n^+" (21) 

is solution of Eq. f|T9|) once setting the dependence a — > a{r). Using once again Eqs. ©-([H]), 
we obtain air) = p{r)/p{r) and hence the most general solution for the problem is given by 

1 — 3- 1 + - / \ 

p = p{m{r),n{r),r) = m pU p, (22) 

where, for simplicity, we have omitted the explicit dependence of m, n, p and p on r, and 
we have adjusted the arbitrary constant n in order to have p = mn in the nonrelativistic 
limit, and p = n*/^ in the ultra-relativistic limit. 

For our purpose, a simple and realistic model for a symmetrically spherical stellar object 
arises from the assumption that the fluid is incompressible: the density po is a constant out 
to the surface of the star, r = R, after which it vanishes. Specifying p{r) is equivalent to 
employ the equation of state, since p{r) can be determined from the hydrostatic equilibrium. 



i. e. from the TOV equations [26|, 



Mr) ( , , ^, (M(r)+47rr^p(r)) 

^ = -(PM +Vir)) , (23) 

and 

^ = 47rrVM, (24) 

for which we have set Newton's constant, G = 1, and M is the total mass of the stellar 
object. 

Integrating Eq. ( l23l) yields 

p{r) [I - 2{M/R){r'/W)]'''' - [1 - 2{M/R)\''' 



Pir) 3 [1 - 2(M/i?)]'/' - [1 - 2{M/R){r^/R^)Y^^' 
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(25) 



Notice that the pressure increases near the core of the star, as expected. Figs. [2ll3] illustrate 
the properties of the relevant thermodynamic quantities for the case where their behavior 
are constrained by Eq. fl22p in the limiting case where the mass of the particle is given by 
its observable value m{R) = tjIq. We compute the ratio between pressure and energy density 
p/p, the density of particles n (dimensionally normalized by ml), and the ratio {p + p)/n 
that characterizes a pseudo chemical potential. 

Indeed, for a star of fixed radius R, the central pressure tends to infinity if the mass 
exceeds 4i?/9, and exceeds the radiation pressure (p = p/3) if the mass exceeds 5-R/18 
(solid hne in the Fig. [2]). 

If we try to squeeze a mass greater than the Buchdahl's limit 4R/9 within a radius R, 



general relativity admits no static solutions 



; a star shrinks to such size till it eventually 



turns into a black hole. This scenario can be modified for astrophysical bodies composed 
by particles subject to the mass varying mechanism. Depending on the analytical mass 
dependence on the cosmological background scalar field modified by the spacetime curvature, 
the stability of equilibrium configurations can be drastically modified, as we shall discuss in 
the following section. 



IV. EQUILIBRIUM AND MASS DEFECT 

Among various types of equilibrium in stellar and cosmological scenarios, the simplest 
one is that in which all processes of reciprocal transformation (creation/annihilation) take 
place at a time scale much faster than the one associated to the flux of external particles 
and the rate of change of the relevant thermodynamic quantities. Such a simplified ver- 
sion of thermodynamic equilibrium is consistent with the stellar structure of our approach. 
The conditions for lumping of such an exotic nonbaryonic matter plus scalar field fluid is 
determined by the total binding energy of the system. 

The form of the Eq. fl2^ suggests that we interpret M as the total (mass) energy inside 
a radius R, including the rest mass-energy Mq, the internal energy of motion W and the 
(negative) potential energy of self-gravitation U. In the relativistic domain, we have 

M = Mo + Mb = 4n 1 r^p{r) dr, (26) 

^0 
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where 



Mo = 47r / r''m{r)n{r)>yA{r) dr. 



(27) 







Following Eq. ([8]), M from Eq. fl26l) is reduced to Mq (no binding energy and no motion) if 
p = 0. The difference between the total rest mass-energy and the total energy corresponds 
to the (positive) binding energy Mg which keeps the stellar structure stable, 



Although it is useful to define M as the total energy, the rest mass-energy Mq is a 
fundamental quantity in the stability analysis. In opposition, the internal energy of motion, 
W, and the (negative) potential energy of self-gravitation, U, are not particularly relevant, 
unless in the Newtonian approximation where the relation Mg ^ U + W holds. The binding 
energy is sometimes referred to as mass defect and it corresponds to the energy which is 
released during the formation of a star from an initially rarefied matter: a typical mechanism 
for compact structure formation. When particles are combined into a bound system, an 
energy which equals the mass defect is emitted in the form of photons, relativistic neutrinos, 
or gravitational waves. For this reason, it follows from the physics of this process that the 
condition for a stable static star originating from diffuse matter distribution is expressed as 



In the Fig. H] we show how the mass varying mechanism does modify the equilibrium 
condition based on the mass defect criterium for relativistic stellar objects. 

Two interesting effects can be described through the mass varying mechanism depending 
on whether the particle mass increases or decreases inwards towards the center. From the 
behavior of the mass defect for compact objects composed by particles with their mass 
exponentially increasing inwards to the center, we notice that stable configurations are 
favored relatively to the other analytical dependencies of m(r). The resulting compact 
structure has binding energy that increases to infinite as the ratio M/R approaches to its 
Buchdahl's limit. When these objects shrink to such a size they inevitably keep on shrinking 
and eventually give origin to a black hole. 

In the absence of the mass varying mechanism, the rest mass-energy Mq is written as 
Mq = mN, where is the total number of particles within the radius R, that is 




(28) 



Mg > 0. 




(29) 
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The stability curve for this case is described by the sohd hne depicted in the Fig. IH for 
which the binding energy has an upper hmit, M/R < 4/9, whether one assumes that the 
pressure at the center cannot go to infinity. It is important to reahze that, in case of mass 
varying particles, the analysis of the stability conditions in terms of N has to be rephrased 
in terms of Mg. 

In this context, decreasing the mass inwards towards the center can also give origin 
to stable structures up to certain limiting values for M/R. In opposition to the previous 
result, the equilibrium conditions are achieved for smaller values of M/R, which restricts 
the existence of neutrino lumps with arbitrarily small masses. The point here is that the 
coupling with the background scalar field is crucial in determining the stability conditions. In 
this situation, the role of the scalar field potential on the mass varying mechanism requires 
specific considerations, which not only depend on the cosmological nature of the scalar 
field (quintessence, cosmon, phanton, Chaplygin gas, etc), but also on the mass generation 
mechanism. 

V. STABILITY OF EQUILIBRIUM CONFIGURATIONS 

The condition of hydrostatic equilibrium is equivalent to the condition of an extremum 
of the total energy of a spherically symmetric stellar object for a given number of conserved 
particles and a given specific entropy. In fact, the relativistic theory of radial perturbations 
for nonrotating equilibrium configurations is well understood. There are two approaches to 
radial stability: a dynamical one, based on the equations of motion and energy properties 
of radial perturbations; and a static one, based on mass-radius curves for certain sequences 
of equilibrium configurations. From the dynamical approach we can find different ways for 
studying the normal radial modes of a relativistic stellar model. From the static approach 
we obtain the minimum energy corresponding to stable equilibrium, and a maximum to 
unstable equilibrium. In the static approach the study of stability does not require additional 
computations. It clearly reveals that general relativity catalyzes radial instability in stellar 
models. 

The analysis of stability conditions from the static approach is similar to one when treat- 
ing mass varying particle systems with thermodynamics governed by Eq. ( I22l) . However the 
analysis has to be performed in terms of Mq and not in terms of A^. 
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First, we consider how the mass of an equihbrium configuration changes with the addition 
of particles, brought from infinity, where the energy is M . In other words, we derive the 
analytical correspondence between M and Mq and, in particular, dM/dMo. The change in 
M should not depend on the position in the equilibrium configuration at which the particle 
is added up. 



It can also be shown that for stable equilibrium configurations, dM/dMo < 1 (39|, |40 |. 
This does not necessarily imply that for equilibrium M < Mq for equilibrium, which is 
an independent condition. In fact, the equilibrium solutions with negative binding energy 
Mg are unstable. Moreover, the stability against small perturbations is lost at the point 
of greatest (positive) binding energy, so that the configurations of negative binding energy 
must lie deep into the unstable region. It is important to realize that the NR analysis, which 
is essentially qualitative, allows for solutions of negative binding energy [soj]. The stability 
condition is reduced to dM/ dA^ < m for the case where the particle masses does not depend 
on the radial coordinate, m(r) = m = ttiq. For this reason, the usual general relativity 
analysis of stability is performed in terms of M instead of A^. In the case of mass varying 
particle systems it gives origin to significant deviations from the equilibrium and stability 
scenarios described in terms of the condition M < Mq. 



A. Stability of mass varying particle lumps 

By considering the above discussed static approach for the study of stability, let us 
describe the stability curves for the considered mass varying particle analytical relationships, 
and compare them with the analysis in terms of the particle number N . In Fig. [5] we plot 
the stability curves, i. e. the relationship between the total mass energy, M, and the rest 
mass-energy, Mq, and its correspondence with respect the analysis in terms of the total 
number of particles contained within R. 

In the Fig. [6] we plot the stability conditions, i. e. dM/dMp and dM/dN for the case 
where m is given in units of itlq. 

The stable equilibrium configurations are achieved if dM/dMg < 1, and not if dM/dN < 

1. 
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B. Neutrino Lumps 



Neutrinos coupled to the dark energy, an ultra light scalar field lead, in its simplest 
formulation, to a number of significant phenomenological implications. Under general con- 
ditions, they materialize in a rather model independent way 



12 



The neutrino mass 



is presumably generated from the vacuum expectation value (VEV) of the scalar field 
(j) whose dynamics is governed by the Universe scale factor a, (p = 0(a), which turns the 
neutrino mass into a dynamical quantity. 

In most of the scenarios, the neutrinos remain essentially massless until recent times. 
Their mass eventually increases close to its present value and their interaction with the 



background scalar field almost ceases 
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161 ]. The energy of the scalar field becomes the 



dominant contribution to the energy density of the Universe and the cosmological accel- 
eration ensues. For the coupled neutrino-scalar field fiuid the squared sound speed may 



become negative - a signal of instability 



251]. The natural interpretation of this instability is 



that the Universe becomes inhomogeneous due to nonlinear fiuctuations of neutrinos, which 



eventually collapse into lumps 



4l| 



The question here is if the fraction corresponding to the neutrino energy contribution 
is due to an isotropic and homogenous NR {p ~ 0) distribution of particles, i. e. Mq, or 
if it is due to a gas of weakly interacting particles {p > 0) which lump in an expanding 
Universe, that is M. The rates given by M^/M for different scenarios of mass varying 
mechanism leads to a reasonable estimative for the corrections on the absolute values of the 
cosmological neutrino masses. 

So far, the known cosmological neutrino mass predictions are performed for an isotropic 



and homogenous NR distribution of particles 



42 



43 



44 



45| , which in the absence of internal 



interactions and gravitational forces yield a total energy (per equivalent lump) equal to Mq. 
It would be interesting to quantify the contribution of neutrino lumps in the scenarios of 
cosmological neutrino masses. 

For the case where the total neutrino number is assumed to be conserved (clearly 
an approximation), one can define the apparent (measured) value for neutrino masses as 
(/i) = M/N and the realistic (expected) corresponding value as (m) = M^/N so that 

M - ^ (301 
Astrophysical objects composed by neutrinos which similar features have also been studied 
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in Refs. [32, |46|, |47|. 

Let us then consider that the energy density of the Universe involves weakly interacting 
neutrinos described by the previously discussed static, spherically symmetric solution of the 
Einstein equation with constant energy density. The neutrino mass depends on the value 
of a slowly varying scalar field, m{(f>{r)) ~ m(r), which is responsible for the attractive force 
responsible for the lumping. Independently of the analytical expression of the mass varying 
mechanism, we assume that the neutrino mass is equal to rrio at R, that is on the surface of 
the spherical lump where p = 0. 

In the Fig. [7] we represent the rate (m) / (/i) for different mass varying scenarios aiming 
to obtain realistic values for the neutrino masses. We observe a clear correspondence with 
the stability curves of the Fig. m 

Naturally, the smaller is the spherical radius, the more isotropic and homogeneous is the 
cosmological neutrino distribution. For lumps in which the neutrino mass decreases from 
the boundary to the center, we have a lower limit for the rate (m) / (/i) , which coincides 
with the upper limit for the radius of stable equilibrium configurations. For the other cases, 
that is, for static or increasing neutrino masses from the surface to the interior, we get 
stable equilibrium configurations for which the rate (m) / (fi) decreases till vanishing at R. 
From the analysis of equilibrium and stability, we notice that an unlimited quantity of mass 
energy M can be squeezed into a lump when the particle mass increases inwards towards the 
center of the star, which is a preliminary condition to form a black-hole late configuration. 
A crude interpretation of this is that the apparent neutrino mass parameter, (m), assumes 
smaller values in the Schwarzchild exterior solution. Since the analytical dependence of 
the scalar field on r is suppressed, the value of (m) can be interpreted as the cosmological 
neutrino mass. Although somewhat paradoxical, this is a natural consequence of the mass 
varying mechanism. Actually, for these lumps, most of the mass is concentrated at the 
center. Assuming that predictions for the neutrino energy density and for the corresponding 
particle number density are unaltered by the lumping scenarios, the neutrino mass value 
vanishes at the boundary (r = R) of this objects if the ratio M/R reaches the Buchdahl's 
limit. The neutrino masses in this cases are expected to vanish. 

On general terms, all stable scenarios of compact lumps lead to a shift towards smaller 
values for the neutrino masses. Comparing two possible scenarios for cosmological neutrinos 
subject to a coupling to the background scalar field, one with an isotropic and homogenous 
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energy distribution, vis a vis another whose perturbations result in stable neutrino lumps, 
the latter correspond to lower masses at present, that is, < 0.07 eV. 

VI. CONCLUSIONS 

In this work we have shown that the presence of an interaction between nonbaryonic 
matter, either neutrinos or dark matter, to a cosmological background scalar field associated 
to dark energy, has interesting implications, besides the changes on the evolution of energy 
density components of the Universe. 

In our approach, we have assumed that any kind of dynamical mass behaves like an 
additional extensive thermodynamic degree of freedom, depending on the background scalar 
field. Under adiabatic conditions, the thermodynamic pressure and its explicit dependence 
on the curved space radial coordinate could be expressed as p(r), and to be computed from 
the TOV equations for the hydrostatic equilibrium 



27] . This justifies the assumption 
of an explicit analytical dependence of the variable mass, m{r). Consequently, a connection 
among the energy density the particle variable mass and the particle density is established. 
Furthermore, the analysis provides the necessary conditions to examine the equilibrium and 
stability of static, spherically symmetric compact structures 

For cosmological scenarios which admit some kind of mass varying mechanism, the insta- 
bilities associated with an imaginary speed of sound in the NR regime do play a role in the 
structure formation. Effectively, the scalar field mediates an attractive force between the 
nonbaryonic particles and may lead to the formation of matter lumps. Actually, this would 
induce the combined fluid to form compact structures which would behave like cold dark 
matter. The detection of these compact objects could be made through their gravitational 
potential, or through their effect on baryons. Therefore, the existence of nonbaryonic matter 
umps coupled with baryonic structures like stars or cold matter, should not be discarded 



48|. 



We have also established a direct connection between the size of neutrino lumps and the 
neutrino masses. 

One interesting property of our solutions is that small modiflcations on the potential of the 
scalar fleld does not affect signiflcantly the analysis of equilibrium and stability conditions. 
Our results reinforce the argument that the Higgs [49| and the neutrino ^] sectors are 
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possibly the only ones where one can couple a new SM singlet without upsetting the known 
phenomenology. We believe that our proposal is a further concrete step in this respect. 
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FIG. 1: Evolution for the variable mass dependence/W^he radius r. Its dependence is supposed to 
be dictated by a spherically symmetric distribution of matter due to the coupling with a background 
scalar field (j){r). 
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FIG. 2: Thermodynamic variables for a system composed by mass varying particles, where the mass 
is assumed to be an extensive thermodynamic degree of freedom. The plots depict the case for 
which the mass of the particle is given by its stationary observable value m{R) = ruo on the surface 
boundary. Notice that increasing gray level corresponds to increasing values of the thermodynamic 
variables, for which the boundary values have been marked for M/R = 5/18, the "soft" Buchdahl's 
limit for which p = 3p at R = 0. 
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FIG. 3: The ratio between pressure and energy densl't^'' p/ p, the density of particles n (dimensionahy 

normahzed by ml) and the ratio {p + p)/n. We assume the "soft" Buchdahl's hmit for which p = 3p 

at i? = 0. This corresponds to sections of graphs marked on Fig. [2l 
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FIG. 4: Mass defect (binding energy) for several mas^^^ying particle functions for compact objects 
with uniform energy density. Given that depends on the radial coordinate r for generic classes 
of curved spaces, we quantify the modifications due to our approach by some mass dependencies 
on r as illustrated in the box, where y = r/R. 
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FIG. 5: Stability curves for compact objects in correspondence with the mass varying mechanism 
prescription. We plot M/R x Mq/R - the right analysis element - and M/R x N/R - the wrong 
analysis element - for the case where m is given in units of mo- They are all in correspondence 
with the graphs of Fig. [6l 
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FIG. 6: Stability condition for compact objects in correspondence with the mass varying mechanism 
prescription. We plot dM/dM^ - the right analysis element - and dM/dN - the wrong analysis 
element - for the case where m is given in units of mo- Curves are all in correspondence with the 
graphs of Fig. [3 
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FIG. 7: Modified predictions for the absolute va neutrino masses due to the lumping of 

neutrinos. For all stable scenarios the absolute neutrino mass predictions are shifted towards 
smaller values. 
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